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We have used Monte Carlo simulation to calculate the 
shot noise intensity Si(ui) at 2D hopping using two models: a 
slanted lattice of localized sites with equal energies and a set of 
localized sites with random positions and energies. For wide 
samples we have found a similar dependence of the Fano factor 
F = Si(0i)/2el on the sample length L: F oc L~ a where a = 
0.85±0.02 and a — 0.85±0.07 in uniform and random models, 
respectively. Moreover, at least for the uniform model, all the 
data for F as the function of sample length L and width W 
may be presented via a single function of the ratio W/L 13 , 
with (3 — 2a — 1 ~ 0.7. This relation has been interpreted 
using a simple scaling theory. 



Shot noise at electron transport has been the subject of 
intensive experimental and theoretical research lately (for 
a recent review see, e.g., Ref. 0), because it may provide 
important information about nonequilibrium properties 
of conductors, unavailable from other transport charac- 
teristics. Another motivation for studies of shot noise is 
its direct relation to electric charge discreteness. Namely, 
the smallness of the spectral density of current fluctua- 
tions at low frequency, S7(0), in comparison with the 
Schottky value of 2el, is a necessary condition for quasi- 
continuous charge transfer |^||. Such "sub-electron" 
transfer through conductors with sufficiently high resis- 
tance R and low stray capacitance C may make possible 
several resistively-coupled single-electron devices insensi- 
tive to background charge randomness Q. In this con- 
text, hopping conductors are very promising, so that the 
development of understanding of shot noise in such con- 
ductors seems to be an important task. 

However, though the basic theory of hopping conduc- 
tivity is well developed ||, until recently little had been 
known about noise at hopping. Few publications we were 
aware of had been devoted to narrowband, 1 / /-type noise 
(see, e.g., Ref. || and references therein) rather than 
broadband fluctuations such as shot noise. This is why 
in the recent work of our group M a detailed theoretical 
study of broadband current fluctuations at ID hopping 
was carried out (on the foundation of prior important 
work on statistics of the so-called Asymmetric Simple 
Exclusion Process (ASEP) model ||). 

For uniform, linear ID arrays the low-frequency noise 
depends on the boundary conditions (namely, the fill- 
ing factors /e,/h of the edge sites) and may or may 
not be dominated by boundary bottlenecks. In the for- 



mer case, the Fano factor F — Si(0)/2el tends to a fi- 
nite value of the order of 1 (e. g., for ft = /r = f, 
T = and negligible Coulomb interaction, F = |1 — 2/|) 
i.e., shot noise suppression is insignificant. In the ab- 
sence of boundary bottlenecks (e. g., if /l = frt. = 1/2), 
the Fano factor tends to zero at large number of hops 
N, but only as N~ x / 2 , i.e. much slower than in ID ar- 
rays of tunnel junctions where F — l/N far enough 
from the Coulomb blockade threshold |J|. (This behav- 
ior has been explained [Q using a simple scaling theory 
which also explains other features, i.e. the frequency de- 
pendence St(lj) cx lj^ 1 / 3, in an intermediate frequency 
range.) Nonuniformity of ID hopping systems decreases 
the noise suppression, bringing the Fano factor closer to 
the Schottky value F = 1. 

The goal of this paper is to show that the ability of 
electrons to circumvent transport bottlenecks at 2D hop- 
ping leads to a qualitatively different situation. Namely, 
in sufficiently long and broad samples the shot noise may 
be suppressed quite considerably: 



F(xL~ 



(1) 



where L is the conductor length and a « 0.85, even in 
ultimately nonuniform conductors. 

We have employed the usual Monte-Carlo simulation 
technique (see, e.g., Ref. 0) to analyze two different 
models, so far both without Coulomb interaction and at 
vanishing temperature: 

- Model A: hopping between sites with random local- 
ization energies, randomly distributed over a 2D sample, 
and 

- Model B: hopping on a uniform, slanted lattice with- 
out site energy fluctuations. 

In both models, each site may be occupied by just one 
electron, and the rate of (inelastic) transitions between 
the sites is described by the usual formula (see, e.g., Eq. 
(4.2.17) in Ref. @): 



A 



1 - exp(-eij/T) ' 



(2) 



corresponding to the constant density of phonon states. 
Here is the electron energy gain during the hop i — ► j; 
in the absence of Coulomb interaction between electrons 
this gain can be expressed as 



(ej - ej) - eE(x l - Xj), 



(3) 



where E is an external electric field applied along the x 
axis. For relatively short samples special care should be 
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FIG. 1. The nonlinear dc conductivity a — I/WE (nor- 
malized by (Jo = Io/WEo), as a function of electric field E 
for Model A. The straight line is the analytical fit discussed 
in the text. Inset: the same data on a linear scale. 

taken to adequately describe electron transfer between 
the electrodes and the edge localized sites. After exper- 
imenting with various options, we have concluded that 
the same expression (||) may be used to describe this 
transfer, without creating unphysical bottlenecks at the 
electrode-sample interfaces ||. 

In our main Model A, single-particle site energies are 
distributed randomly within a broad energy band, with a 
constant 2D density of states D, and site positions Xi,yi 
are randomly distributed within a rectangular sample of 
length L and width W . The rate amplitude A is an 
exponential function of the intersite distance : 

A = A exp(-rij/a), (4) 

where a is half of the localization radius. All the re- 
sults have been averaged not only over a sufficiently 
long time period, but also over a set of random sam- 
ples with the same global dimensionless parameters L/a, 
W/a, and eEDa 3 (parameter Aq just determines the 
scale Iq — eAoW/ Da 3 of the total current). Such averag- 
ing requires considerable computer resources; the calcu- 
lations have been performed on IBM SP parallel super- 
computer. 

Figure 1 shows the numerically calculated nonlinear 
dc conductivity a = I /WE as a function of electric 
field E for sufficiently long and wide samples. Depend- 
ing on E, we have simulated samples of area L x W 
ranging from 120a x 30a up to 300a x 120a to keep 
the number of "active" sites N s ~ 1500 approximately 
constant (the growing error bars at lower E are due 
to larger fluctuations from sample to sample). The ab- 
sence of significant dependence of a on the sample size 
was being checked. For low electric fields E < 0.02£ , 
(E = 1/eDa 3 ) the current follows closely the depen- 
dence I/E (x exp[-C(i?o / E) 1 / 3 ] expected for 2D variable 
range hopping |5pC[| in the activationless ("high field") 
regime. The best fit (straight line in Fig. 1) gives the 
numerical constant C « 1.27. (This number is to be 
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FIG. 2. Dependence of the Fano factor F in Model A (av- 
eraged over 32 sample realizations) on the sample length L for 
E — 0.035-Eo. Circles: W = 5a, asterisks: W = 10a, squares: 
W = 20a, diamonds: W — » oo. Error bars show the standard 
deviation of the mean. 

compared with the value C = 1.02 following from an- 
alytical calculations in Ref. The minor deviation 
from the analytical dependence at E > 0.02i?o is possibly 
due to multiple, well-branched percolation paths which 
are not considered in the usual theoretical treatment of 
high field hopping. 

Figure 2 shows the Fano factor (averaged over 32 
sample realizations) as a function of the sample length 
for E = 0.035-Eo and for tree values of sample width: 
W = 5a (circles), W = 10a (asterisks) and W — 20a 
(squares). At this field the average hop length is r « 3.3a 
(with r.m.s. projections r x rj 3.1a, r y w 1.9a) so that the 
factor exp(— r/a) is still small. Most importantly, we 
see that shot noise is suppressed considerably (F -C 1) 
in wide and long samples. As a function of the sample 
length L at fixed width W, the average Fano factor first 
decreases following Eq. (0) (solid line in Fig. 2) and then 
at certain length deviates up from this dependence. The 
deviation starts at larger L for wider samples. For nar- 
rower samples one can see the saturation of the average 
Fano factor at large L; simultaneously the width of the 
Fano factor distribution grows significantly. The average 
Fano factor also decreases with the sample width at fixed 
length, obviously saturating at large W since the electron 
transport in remote parts of a very wide sample is un- 
corrected. Having performed the calculations at fixed 
length for several widths (not all the results are shown in 
Fig. 2) we have extrapolated the Fano factor dependence 
to infinitely wide 2D samples. As a function of L, these 
results (shown in Fig. 2 by diamonds) closely follow Eq. 
(0) with a = 0.85 ± 0.07. The power-law dependence 
F(L) has been observed within one order of magnitude 
range of L. The longer samples have not been studied due 
to computer limitations (as an example, calculations of 
a single point L = 160a, W = 20a in Fig. 2 has required 
380 hours of total CPU time). 

In order to verify the shot noise suppression for larger 
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FIG. 3. Fano factor dependence on the length N for hop- 
ping on a uniform slanted square lattice (Model B, inset) with 
width M equal to 2, 5, and N. 

set of lengths and widths, we have used the simplified 
Model B in which (N — 1) x M localized sites with equal 
energies are arranged on a uniform slanted square lattice 
(see inset in Fig. 3) p2[ . In accordance with Eq. (|2|), at 
T = the transport is unidirectional, and transfer rates T 
between all the internal neighboring sites are equal. For 
the links from the left electrode to the nearest internal 
sites we have selected the rates 2r/^ (the factor of 2 
reflects two "channels" per internal site) while the rates 
of hopping onto the right electrode are 2r(l— fa). For the 
numerical analysis we have chosen the case /l = fn = 0.5 
in which, similarly to ID model, there are no boundary 
bottlenecks. Since Model B does not require averaging 
over different random realizations, it may be studied with 
much better accuracy using the same computer resources. 

Figure 3 shows the Fano factor as a function of the 
array length N for several values of width M . Again, 
we see a strong shot noise suppression. For sufficiently 
wide samples the suppression follows Eq. (|l]) (where now 
L should be replaced by N), with similar exponent as in 
the wide random samples: a = 0.85 ± 0.02. On the other 
hand, for fixed M and sufficiently long N the suppres- 
sion power approaches 0.5, i.e. the same value as for ID 
hopping §. 

There is numerical evidence (see Fig. 4), as well as 
scaling arguments (see below), that for N ^ 1 and M 3> 
1 the crossover between these two asymptotic laws may 
be parameterized in the following way: 

F(N,M) = N-^M-ig(N^/M), (5) 

where (3 — 0.70 ±0.04, and the function g(x) is shown in 
Fig. 4: 

. f const, x > 1 
g{x) K ( x-V 2 , x « 1. (6) 

Checking if the data from Fig. 2 for the random Model 
A may also be collapsed on the similar universal curve 
we have found a reasonable fit for relatively wide sam- 
ples (see inset in Fig. 4) for the following replacements: 
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FIG. 4. The data from Fig. 3 (circles: M = 2, triangles: 
M = 5, and stars: M — N) collapsed onto universal curve 
g(x) using parameterization x = N ' 7 /M and g = (NM) 1/2 F. 
Lines connecting points are just guides for the eye. The inset 
shows the same data with added points for W = 20a (squares) 
from Fig. 2 using the rescaling N = L/5.9o and M — W/lla. 

N = L/5.9a, M = W/lla (for the particular field 
E = 0m5E o ). 

We will start the interpretation of our findings from 
the uniform Model B. Similarly to the ID ASEP model 
, in the absence of lateral boundary effects due to finite 
M, and at = = / the probability of any charge 
configuration is expected to be the same as if each site 
had independent occupation with probability p = f [fl3[ . 
As a result, dc current between any two neighboring sites 
should equal Tp(l — p), so the total dc current is 

I = 2MTp(l - p). (7) 

Following the arguments of Ref. we obtain only a 
minor suppression of low- frequency shot noise, F w |1 — 
2p|, in the case p ^ 0.5 (the noise significantly decreases 
at frequencies u> > oji ~ r|l — 2p\/N). However, when 
the coupling with electrodes is strong enough, /l > 0.5, 
Jr < 0.5, the half-filling is expected inside the array, 
p = 0.5, and the Fano factor can indefinitely decrease 
with the array length N. In this case, for sufficiently 
large N (narrow array) we may repeat all ID scaling 
arguments of Ref. fjj based on Eq. (0) , and arrive at the 
following estimates: 

F - (NM)- 1/2 , lui - TN~ 3/2 M- 1/2 , (8) 

for the Fano factor and the saturation frequency u>i, 
above which the dependence Si(ui)/2el ~ (w/r) -1 / 3 ^ -1 
is expected. 

Obviously, this result should eventually fail if we start 
increasing the width M for fixed N, because F cannot de- 
pend on M in the limit of wide array (since the transport 
in remote parts of the array is uncorrelated and so both 
Sj and I are additive). Denoting the crossover width as 
Mq we get the estimate F ~ (NAIq)^ 1 ^ 2 for wide ar- 
rays. It would be natural to expect M cx N, however, 
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the numerical results (Fig. 3) indicate the power-law de- 
pendence, Mq oc N@ , with (3 being a phenomenological 
parameter. If this assumption is true, we obtain the fol- 
lowing estimates: 

F~iv-( 1+ «/ 2 , ~ riv-( 3 +«/ 2 , (9) 

for wide arrays, M 3> N 13 . Our numerical result, a = 
0.85±0.02, then leads to the value f3 = 2a- 1 = 0.7±0.04. 
For intermediate widths it is natural to suggest that the 
crossover is governed by some function of the ratio x = 
N° /M alone. Thus we recover the behavior described by 
Eqs. dM|) and illustrated by Fig. 4. [Actually, at this 
stage we cannot rule out the possibility that the function 
g in Eq. (||) also has a weak dependence on M that would 
lead to either smoothing or sharpening of the curve in 
Fig. 4, leaving however the asymptotes (|J) intact.] 

If we apply the same scaling arguments which have led 
to Eq. (||), to smaller blocks of size x N%, we obtain 
the noise frequency dependence 

S J (w)/2eJ~o;-TJV- 1 , 7 = (1 - 0)/(3 + 0), (10) 

for wide arrays at intermediate frequencies, toi <C u) <C L. 

This power-law dependence, Si(ui) oc cj~ 08 , was con- 
firmed numerically with the accuracy of the exponential 
factor about ±0.007. The same frequency dependence 
should be expected for the narrow arrays, M <C at 
frequencies higher than ~ FM~( 3+ ' 3 )/ 2 ' 3 , while at lower 
frequencies Si(u>) oc w -1 / 3 , as discussed above. 

Now, let us turn back to our main, random Model A. 
According to the percolation picture of hopping H , the 
conductivity of a sample is determined by a "percolation 
cluster" , essentially a network of sites connected by the 
most probable hop paths. The percolation cluster may 
be divided into blocks of a certain size, such that the 
aggregate characteristics of each block (e.g. the average 
current) are nearly equal, even though inside each block 
the sample is highly (exponentially) nonuniform. Aver- 
age block size in the transport direction for the random 
model A may be determined by mapping 2D hopping in 
disordered wide samples onto the uniform model B, for 
M > N, where the shot noise suppression was found to 
be similar. If this interpretation is valid, then the ratio 
L/N = 5.9a obtained from the mapping should be com- 
parable to the correlation length L c ~ f(f/a)" ~ 3.3a 
]14j |. We believe these numbers are reasonably consis- 
tent. 

On the other hand, the model A behavior in narrow 
samples (W < 10a) is quite different from that in the 
uniform model: instead of going down with growing TV, 
the Fano factor saturates. Simultaneously, the statistics 
of F becomes much wider. This behavior is very natural, 
since if the sample is narrower than the block size (for 
our value of electric field, about 10a), the exponentially 
broad distribution of hopping paths within the block is 



revealed and is mapped onto the properties of the sample 
as a whole. 

Our results for wide samples are in good agreement 
with data from a recent experiment [ pTsj in which shot 
noise at hopping was measured in p— type SiGe quantum 
wells. Actually, the experimental F V curve significantly 
differs from that in our Model A. However, our result 
(|l|) for F seem more general. For wide samples with two 
different lengths L\ = 2/zm and Li = 5/zm the Fano fac- 
tor was measured |ll| to equal F\ — 0.43 and F2 = 0.2, 
respectively. This corresponds to Eq. ([[]) with a — 0.84, 
the value which is virtually equal to our result a = 0.85. 
Such perfect agreement is possibly just a coincidence, 
since so far only two experimental points are available. 
Evidently, it would be valuable to have more experimen- 
tal data, in order to verify the shot noise suppression 
power a = 0.85. 

In conclusion, we have numerically investigated shot 
noise suppression at 2D hopping in random samples as 
well as in uniform arrays. Very similar shot noise sup- 
pression has been found for both models for wide and 
long samples. At 2D hopping the Fano factor decreases 
with the sample length as F oc 1/L - 85 . This suggests 
that shot noise suppression at 2D hopping is insensitive 
to details of the hopping process, e.g., the energy depen- 
dence of the hopping rate. If this surmise is true, the 
Fano factor should not be a very strong function of tem- 
perature. It may, however, be substantially altered by 
Coulomb interaction of hopping electrons, as in the ID 
case . Our next plans are to explore the effects of both 
these factors. 
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